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Abstract. The purpose of this contribution is to review some aspects of the loop space 
formulation of pure gauge theories having the connection defined over a Lie algebra. 
The emphasis is focused on the discussion of the Mandelstam identities, which provide 
the basic constraints upon both the classical and the quantum degrees of freedom of 
the theory. In the case where the connection is extended to be valued on a super Lie 
algebra, some new results are presented which can be considered as first steps towards 
the construction of the Mandelstam identities in this situation, which encompasses 
such interesting cases as supergravity in 3 + 1 dimensions together with 2+1 super 
Chern-Simons theories, for example. Also, these ideas could be useful in the loop space 
^*| . formulation of fully supersymmetric theories. 



^ 1. INTRODUCTION 



Gauge theories provide a successful framework to unify the strong, weak and 
electromagnetic interactions in nature. The recent introduction of the Ashtekar 
variables in general relativity has permitted to reformulate the gravitational inter- 
action in the form of a standard gauge theory, plus additional constraints. Gauge 
theories where no matter is included, are fully described in terms of a single field: 
the connection A^, which is the non-abelian generalization of the standard electro- 
magnetic potentials. In general, gauge theories are characterized by being invariant 
under transformations generated by local symmetries, i.e. transformations of the 
connection which may be differently chosen at each point in space time. This free- 
dom has the consequence of introducing arbitrary functions in the dynamics of the 
connection. In this way, only those specific fields constructed from the connection 
which are independent of these arbitrary functions, will be of physical relevance: 



they are called gauge invariant objects. From this point of view one could say that 
taking as the basic field in a gauge theory is an unnecessary complication. This 
observation has prompted many proposals to describe gauge theories starting ab 
initio with only gauge invariant objects. 

In this contribution we will consider the so called loop representation approach, 
to be described with more detail in the sequel. This method has been very suc- 
cessful in the description of non-perturbative aspects of gauge theories and also in 
providing spectacular advancements in the problem of quantizing Einstein gravity. 
The basic underlying idea of the loop representation is that the Hilbert space of 
gauge invariant states can be spanned by states which are labeled by loops. A main 
feature of the loop approach is the redundancy of the gauge invariant degrees of 
freedom: they constitute an overcomplete set which is restricted by the so called 
Mandelstam identities. These identities must be taken into account, either by im- 
posing them over the physical states in a manner similar to the standard Dirac 
constraints, or by explicitly solving them in order to identify the corresponding 
reduced phase space. Even the preliminary problem of identifying the complete set 
of Mandelstam identities for a given physical situation is not completely solved in 
the general case for gauge theories over a Lie algebra. 

In this work we will discuss only theories described by a connection field, i.e. 
those where matter fields are absent. These include pure Yang Mills theories, pure 
2 + 1 dimensional Chern Simons theories and also gravity and supergravity in the 
Ashtekar formulation, for example. From the Hamiltonian point of view, these 
theories are basically on the same footing. Each one is characterized by a specific 
canonical Hamiltonian together with a given set of first-class constraints. We will be 
mostly concerned with the elucidation of the problem of the Mandelstam identities, 
both in the standard case and also when the connection of the gauge theory is 
defined over a super Lie algebra. Besides its own interest, the latter situation 
is also relevant as a first step in the application of loop representation methods 
to fully supersymmetric gauge theories, which involve matter fields as well. The 
discussion will be focused at the classical level, without considering in detail the 
construction of quantum states as functionals of loops. These notes do not have the 
pretense whatsoever of being a review on the subject and from the very beginning 
we apologize to those authors whose work we have involuntarily not cited. 

The general organization of this paper is as follows: in Section 2 we provide 
a brief review of the standard formulation of Yang Mills theories in terms of the 
connection. We recall the Lagrangian formulation as well as the Hamiltonian ap- 
proach, in which the loop representation is based. Section 3 contains the basic 
material necessary to extend the connection to a super Lie algebra object and 
includes the definition of a supermatrix together with those of the fundamental 
invariants under similarity transformations: the supertrace and the superdetermi- 
nant. An introduction to the loop space formulation of gauge theories is presented 
in Section 4, where the basic concepts are reviewed to end up with the introduction 
of the gauge invariant degrees of freedom: the Wilson loops. Also, the construc- 
tion of the wave functionals of the system in the loop space is briefly sketched. 



Section 5 is devoted to the discussion of the overcompleteness of the loop space 
variables and the Mandelstam identities are introduced to account for this redun- 
dancy, both in the classical and quantum case. The so called generic Mandelstam 
identities are subsequently derived from the identities arising from the application 
of the Cayley-Hamilton theorem in the specific representation where the connection 
lives. Section 6 deals with the extension of the previous ideas to the case where 
the connection is defined over a graded Lie algebra. To begin with, we present an 
extension of the Cayley-Hamilton theorem to the case of supermatrices. Besides its 
own interest, our motivation to deal with this mathematical problem is the claim 
that the algorithm developed for the standard case in Section 5 can be extended to 
this situation, thus providing the corresponding Mandelstam identities. The exten- 
sion of the Cayley-Hamilton theorem proceeds in two steps: (i) the identification 
and definition of a characteristic polynomial for supermatrices and (ii) the proof 
that each supermatrix annihilates the null polynomial previously defined. Also, 
examples of characteristic polynomials for supermatrices are given in some simple 
cases. Finally, examples of Mandelstam identities constructed for some simple su- 
permatrices, according to the previous ideas, are presented. We close with Section 
7, which contains some conclusions and a brief outlook of some remaining open 
problems related to the topics discussed in this paper. 



2. THE CONNECTION APPROACH TO YANG MILLS 

THEORIES 

2.1. The Lagrangian Formulation 

Let us consider the Yang Mills theory corresponding to the compact group G, 
which is characterized by its Lie algebra A having the antihermitian generators 
T A , A = 1, . . . ,ra. They satisfy the commutation relations [t a , T b ] = f AB c T c , 

where f AB c are the structure constants of the group. A standard representation 
of these generators is in terms of ordinary matrices, with additional restrictions 
appropriate to the group under consideration. 

The basic object in this formulation is the connection A^(x), which is a covariant 
vector field that lives in the Lie algebra of G, i.e. A^{x) = A^c(,x)T c ', where 
H — 0, 1, 2, 3 is a world subindex. The infinitesimal local transformations SQ(x) 
generated by the group can be written as SQ(x) = 5Q(x)c T c , where 5Q(x)c are 
arbitrary numerical functions. Under such rotation, the connection transforms as 

SqA^x) = d„6e(x) + [A^x), 8e(x)} := V^Qix), (1) 

which generalizes the gauge transformation of electrodynamics. Here we have in- 
troduced also the covariant derivative of any object Mc with one subindex in the 
Lie algebra, in the form 

V li M = d li M + [A li (x), M] , (2) 



where M = M c T c . In components, the above expression reads 

(V I1 M) C = d,M c + f BD c A^ B M D . (3) 

The object which is covariant under the local group of transformations is the 
field strength F pv = F pvC T c defined by 

F^ = d p A u - d u A^ + [An, A v \ , (4) 

which generalizes the electromagnetic field of electrodynamics. Under gauge ro- 
tations, F^ u transforms covariantly , i.e. S@F pu = [F^, 50]. Besides, the field 
strength satisfies the Bianchi identity 

V^F up + V v F p(i + V p F^ = 0. (5) 

Let us emphasize the basic property that the covariant derivative (2) of any covari- 
ant object is itself a covariant object. 

The dynamics of the pure Yang Mills is described by the action 

S = Jd 4 x £ = Jd A x ^Tr(F^F^), (6) 

where C is the Lagrangian density and Tr(T A T B ) = h AB is an invariant symmet- 
rical tensor in the corresponding representation of the Lie algebra. For a semisim- 
ple group (i.e. one having no U(l) invariant subgroups) it is possible to choose 
h AB = —5 AB . The world indices are raised and lowered by the Minskowsky met- 
ric 7]^ = diag(l, — 1, — 1, — 1, — 1). The action (6) is invariant under the gauge 
transformations (1). Finally, the resulting equations of motion are 

£> M f" = o. (7) 



2.2 The Hamiltonian Formulation 

This formulation corresponds to the (1 + 3) splitting of Minkowsky space-time, 
where we select three-dimensional hypersurfaces of constant time. It is convenient 
to introduce the chromoelectric, E i: and chromomagnetic, B\ fields in the following 

way 

Ei = F 0l = Ai- ViA , B { = l -e^ k F jk . (8) 

Here, i = 1,2,3 labels the components of a vector living on the three-dimensional 
hypersurface and the dot denotes the time derivative evaluated on that hypersur- 
face. The covariant derivative T>i corresponds to the definition (2) restricted again 
to the hypersurface, through the conecction A^ Introducing the canonical momenta 



and denning the Hamiltonian density H = U l cAic — C, we obtain the action 



S = Jdtjd 3 x (n'cAc - \{E 2 + B 2 ) + Aoc(A^)c) , (10) 

where E 2 = E l c E % c > and analogously for B 2 . The above expression arises 
after an additional integration by parts, which can be performed in virtue of the 
following properties of the covariant derivative: (i) T>i(MN) = {T>iM)N +M{T>iN) 
and (ii) Tr(T>iM) = diTr(M). The action (10) tells us that the phase space 
variables are the coordinates Aic(x, t), together with their conjugated momenta 
W D (y,t) = Ej D (y,t), which satisfy the standard Poisson brackets 

{A iC (x,t), W D (y,t)} pB = 5 3 (x-y)5^5 CD , (11) 

at equal times. Also, we infer from (10) that A oc (x) are Lagrange multipliers 
leading to the Gauss law constraints Gc = iT^iE l ) c ~ 0, where Gc is such that 
g = ViE* = Q c T c . In particular 

g c = d^ c + f BD c A lB K D . (12) 

Applying the standard Dirac method for constrained systems [1,2], we can verify 
that there are no secondary constraints. Also, one obtains that Qc are first-class 
constraints, which generate gauge transformations on the three-dimensional hyper- 
surface. A compact way to show these properties is by introducing the averaged 
constraints 

Ga = J d 3 x A c g c , (13) 

where A c \x) are arbitrary functions of position. In this way, the Poisson brackets 
of the constraints turn out to be 

{Ga, Ge} PB — ^Axe, (14) 

where (A x 0)^ = 2f BC \AbQc- The group indices are lowered and raised with 
the metric h BC together with its inverse hsc- The fact that the constraints Ga 
generate gauge transformations can be seen from the following calculation 

S A A iC (x,t) := {A iC (x,t), Gsa} pb = -(ViSA) c , (15) 

which reproduces the spatial part of the transformation (1). The precise relation 
between the gauge symmetries obtained in the Hamiltonian formulation and those 
appearing in the Lagrangian formulation can be found in Ref. [3]. 



The canonical quantization of the Hamiltonian version of a gauge theory follows 
the standard steps: 

(i) The canonical variables are promoted to the range of hermitian operators: 

A iC — * A iC , E iC — > E iC . 

(ii) The Poisson brackets algebra is turned into a commutator (anticommutator) 
algebra, according to the statistic of the involved fields. In our case the resulting 
commutators are 

A iC (x,t), E jD (y,t) =ih5' i {x-y)5 ij 5cD, (16) 

with all others been zero. 

(iii) A possible way of realizing the above algebra is in the "coordinate " repre- 
sentation, where 

A iC := A tC (x, t), E lC := ik — (17) 

dA iC {x,t) 

The wave function of the system is then a functional of the connection: *tf(A) and 
we can calculate the action of any operator upon it. The basic functional derivative 
is defined by 

(iv) The realization of the quantum operators may lead to ordering problems. 
In particular, the first-class constraints Qc must be realized as operators which 
correctly close under commutation, in order to have a consistent theory. If it not 
possible to do so, we say that an anomaly appears. The first-class constraints are 
subsequently imposed as null operator conditions upon the wave functions repre- 
senting the physical states 

Qc = 0. (19) 
In our case, the operator expression for the Gauss law constraint is 

g e = ihj d 3 x (V t Q) c y-t- (20) 
J 6 A iC 

and the condition (19) simply means that the physical wave functions must be 
gauge invariant functionals of the connection. 

(v) There must be a well defined scalar product which allows for the realization 
of hermitian operators together with the calculation of the transition probabilities 
that permits us to make predictions regarding the observables of the theory. 

(vi) Finally, the dynamics is contained in the wave function ^(A, t), which sat- 
isfies the Schroedinger equation 



(21) 



A corresponding Hamiltonian formulation can be also developed for Chern- 
Simons theories, which are connection theories defined in a space-time having an 
odd number of dimensions. These topological field theories are defined by a La- 
grangian density given by the In — 1 form Qcs defined by dQcs — Tr(F n ). Here, 
F is the 2-form corresponding to F — \F llv dx ll dx v , where F^ v is the field strength 
given in Eq. (4). In the case of 2 + 1 dimensions, the Lagrangian action for the 
Chern-Simons theory is 

S=\Tr [ (dA — ^-A A A) A A, (22) 
2 J m V 3 / 

where the integration is made over a three-dimensional manifold Ai and A = A^dx^ 
is the connection 1-form. 



3. YANG-MILLS THEORIES WITH A CONNECTION 
DEFINED OVER A SUPER LIE ALGEBRA 

Super Lie algebras are a special case of graded Lie algebras, which appear in the 
description of supersymmetry. This new type of symmetry arises when attempting 
to provide a unified description of bosons and fermions. In this case, the basic object 
will be a multiplet incorporating both kind of fields. The allowed "rotations" of 
the theory will mix the components of such multiplet: i.e. bosons and fermions. In 
this 

section we will present only a brief review of some underlying ideas in super- 
symmetry, which are relevant to the construction of supersymmetric connection 
theories [4]. 



3.1. Grassmann Numbers 

Let us start from the "classical" version of supersymmetric theories, which re- 
quires the introduction of Grassmann numbers to represent the algebraic properties 
of fermionic fields. The simplest way to motivate these numbers is by starting from 
the quantum description of two independent fermionic operators / and g, defined 
by the following anticommutators 

{/, / f } = 1, {/, /} = 2/ 2 = 0, {/t, gi} = 0, {/, gt} = 0, (23) 



with analogous expressions obtained interchanging / and g. Consider a realization 
of the above anticommutation relations in terms of coordinates and derivatives in 
a manner analogous to the holomorphic representation for the standard harmonic 



oscillator. To this end, let us introduce two independent coordinates 9 and rj and 
define the following realization of the above operators 

<eri\p = 0<0ri\ < 6 rj\f = < 6 rj\, 

< 9 ri\g ] = ri < ri 9\, < 9 rj\g = ^- < 9 rj\. (24) 

The third relation in Eq.(23) requires 9r)+r)9 = 0, while the second relation demands 
that each of the above numbers must have zero square, i.e. 9 2 = = rf. The 
last relation in Eq.(23) says that j^r\ + r]-^ = 0, which require that the derivative 
operator anticommutes with the coordinates. The above realization of the fermionic 
operators must be supplemented with a scalar product which guarantees that in fact 
0t = -4. For the independent coordinates 9,9*, for example, such scalar product is 
given by 

(*,A) = J d9d9* tf* e 9 * 9 A, 
where the integration over the Grassmann variables is defined by 

J d9 = 0, Jd99 = l, (26) 

and analogously for the complex conjugated variable 9*. 

The "numbers" 9, 9* and rj, r/* that satisfy the above properties are called odd 
Grassmann numbers and provide a "classical" description of fermionic degrees of 
freedom, in the same way as complex numbers allow for the realization of bosonic 
degrees of freedom. The product 9r\ = —r\9 is called an even Grassmann number 
because it commutes either with 9 or rj and also has zero square. Thus, the set of 
all commuting numbers is augmented from the complex numbers to the set of even 
Grassmann numbers. The reader is encouraged to look into Ref. [5] for a detailed 
description of the properties of Grassmann numbers. 

In this way, a unified description of bosonic and fermionic degrees of freedom will 
start from a multiplet, say, 

\u>=[q\ q 2 ,..., q m ;9\9\...,9 n ] T , (27) 

which contains bosonic degrees of freedom q\ i = 1, . . . m, described by even Grass- 
mann numbers (complex numbers in particular), and fermionic degrees of freedom 
9 a , a — 1, . . . , n described by odd Grassmann numbers. The superindex T denotes 
standard transposition, in such a way that the array (27) is a column. 



(25) 



3.2. Supermatrices 

The natural operators acting on the state \u > are linear transformations M that 
produce a "rotated" state \v >= M\u >, which components preserve the even/odd 
character of each entry in the multiplet. Such object is the (m + n) x (m + n) array 



where the corresponding blocks have the following properties: (i) A and D are 
respectively m x m and n x n matrices with all entries been even Grassmann 
numbers, (ii) B and C are respectively m x n and n x m matrices with all entries 
been odd Grassmann numbers. The array M is called a supermatrix and the 
corresponding action upon the state \u > is called a supersymmetry transformation. 

Let us observe that the addition and the multiplication of two supermatrices ( 
with both operations defined in the standard way) is again a supermatrix. Since 
the unit matrix X is also a supermatrix, the inverse of a supermatrix can also be 
defined in the usual manner. 

An important topic in the study of supersymmetry is the construction of the 
supermatrix invariants under similarity transformations, which generalize the basic 
concepts of the trace and the determinant in classical linear algebra. The super- 
trace, denoted by Str, is defined by 

StrM := TrA - TrD, (29) 

where the relative minus sign is required in order to guarantee the cyclic property 
Str(MiM 2 ) = Str(M 2 Mi), under the presence of odd Grassmann numbers in each 
supermatrix. The definition of the superdeterminant, denoted by Sdet, follows the 
same pattern of the classical case and it is given by 

SdetM := exp(Str(lnM)). (30) 

This expression can be written in infinitesimal form as 

5ln(SdetM) = StriM^SM), (31) 

subject to the boundary condition SdetX = 1. The equation (31) is a condensed 
way of writing the partial derivatives of the superdeterminant with respect to the 
supermatrix entries, in terms of the inverse supermatrix M _1 . The superdetermi- 
nant satisfies 

Sdet(MiM 2 ) = Sdet(M 1 )Sdet(M 2 ). Using these properties one can find the 
general expressions 

r,, det(A- BD~ X C) detA 

Sdet ^ = detD = iet[D .cA-^y (32 > 

where we assume the existence of A' 1 and D^ 1 . The determinants appearing in 
Eq. (32) are defined in the usual way because they involve only even (commut- 
ing) Grassmann numbers. The equivalence among the two forms of writing the 
superdeterminant in Eq.(32) is proved in Ref. [6]. 



3.3. The Superconnection 



Motivated by the above discussion, it is possible to construct connection theories 
where is extended to be a (m + n) x (m + n) supermatrix 

A, = A, C T C + ^j, c S c , (33) 

which incorporates bosonic (A^c) and fermionic (Vv) gauge fields that are even 
and odd Grassmann functions of position, respectively. The generators T° and S c 
are purely numerical (complex) matrices which have the generic structure 

t=c 1), s=a t). (34) 



dj' \c 

in the same block form of Eq.(28). They are labeled as even and odd generators, 
respectively. Under an infinitesimal transformation generated by the supermatrix 
5Q = 5QcT c + SQ C S C , the connection transforms in the form given in Eq.(l). In 
order that this transformation be closed in the algebra generated by T° , S c it is 
necessary to require the following (anti) commutation relations among the genera- 
tors 



T , T 



= f AB c T c , [T A , S b ] = g Ab c S c , {S a , S b } = h ab c T c . (35) 

The structure (35) is called a super Lie algebra and it is characterized by the 
presence of anticommutators among the fermionic generators besides the standard 
commutators for the remaining cases. The anticommutators are induced from the 
basic commutator in Eq. (1) due to the anticommuting properties of the fermionic 
components of the multiplet under consideration. We must emphasize that the 
functions which multiply an odd generator in each multiplet , like W or SQ C , 
for example, are fermionic in character , so that they anticommute among each 
other, being realized by odd Grassmann functions of position. The quantities 
f AB CiQ Ab ci h ab c are the numerical structure constants of the super Lie algebra. 

The field strength is now a supermatrix and together with the covariant 
derivative retain their definitions given in Eqs. (4) and (2) respectively. Finally, 
the action for a pure Yang Mills supersymmetric theory is given by an equation 
analogous to (6), 

S = J d 4 x C = J d 4 x ^StriF^F^), (36) 

where the trace has been replaced by the supertrace. The canonical analysis and 
quantization of this theory follows analogous steps to the standard case, except 
that now fermions have been included in the connection. 

This ideas can be extended to the construction of supersymmetric Chern-Simons 
theories. The formulation of gauge theories where the connection is defined over 
a super Lie algebra encompasses some interesting cases like 2 + 1 dimensional 
supergravity with cosmological constant [7] and also the standard supergravity 
in 3 + 1 dimensions, in terms of the Ashtekar variables [8]. The case of global 
supersymmetry is not included in this way. 



4. THE LOOP APPROACH TO GAUGE THEORIES 



The description of gauge theories only in terms of gauge invariant objects, like 
the Wilson loops for example, was initiated by the work of Mandelstam [9]. Other 
earlier references on the subject are [10], [11], [12], [13]. The loop representation, 
which is a quantum Hamiltonian representation of gauge theories in terms of closed 
curves (loops) was subsequently introduced in Refs. [14], [15]. Roughly speaking, 
the dynamical variables in this method are the generalization of the Aharanov- 
Bohm phase of electrodynamics to the general non-abelian case: the Wilson loops. 
In the abelian case, this phase is basically the integral of the connection around a 
closed loop and it is a gauge invariant object. In this way, one considers all possible 
Wilson loops in a three dimensional surface of constant time, as the degrees of free- 
dom of the system. Thus, one shifts from the space of connections to the space of 
loops so that the wave function of the system, originally a functional of the connec- 
tion, turns out to be a functional of loops. Besides allowing for the construction of 
gauge invariant objects, loops provide a natural geometric framework to describe 
gauge theories and gravitation. This approach is a non-local description of the 
dynamics and it is well suited to the discussion of non-perturbative effects. One of 
the central problems in this method is to find and independent set of loop-space 
degrees of freedom, because as we will see in the sequel, the loop-space variables 
form an overcomplete set. 

A very nice example of the properties of this formulation in the case of stan- 
dard gauge theories is given in a paper by Brugmann [16] where the lattice 
gauge theory for the group SU(2) is completely solved within this scheme, and 
compared very favorably with alternative solutions. In the author's words: 
" the eigenvalue problem for the Hamiltonian of SU(2) lattice gauge theory is 
formulated in the loop representation, which is based on the fact that the physical 
Hilbert space can be spanned by states which are labeled by loops. Since the inner 
product between loop states can be calculated analytically, the eigenvalue problem 
is expressed in terms of vector components and matrix elements with respect to the 
loop basis. A small-scale numerical computation in 2 + 1 dimensions yields results 
which agree with results obtained from other methods. " 

Another success of the loop space approach can be found in its application to 
Quantum Gravity [17], [18]. The use of the Ashtekar variables [19], allows to 
rewrite Einstein gravity as a connection theory, in terms of a selfdual connection 
which lives in the Lie algebra of SU(2), plus some reality conditions. The structure 
of the constraints in such formulation is that of the corresponding 

gauge theory plus some extra constraints related to the invariance of gravity un- 
der diffeomorphism in 3 + 1 dimensions. Using the loop space approach, it has been 
possible for the first time to find explicit solutions to the diffeomorphism constraints 
of the complexified theory in terms of functionals of knots. The reality conditions, 
necessary to get back to real Einstein gravity, subsequently provide a definition of 
the required scalar product that leads to an interpretation of the quantized the- 
ory. As we will see in the sequel, loop states automatically solve the local SU{2) 



gauge invariance of the theory. In this way, the three dimensional diffeomorphism 
invariance of Einstein gravity can be realized by labeling the states by knot classes, 
which are just the equivalence classes of loops under diffeomorphism. Finally, there 
remains only one constraint to solve, which is satisfied in terms of superposition of 
knot states. 

In this work we present only a brief introduction to the loop space approach of 
gauge theories. Detailed versions, including applications to quantum gravity, can 
be found in Refs. [20], [21], [22], [23], [24]. 

4.1. Open paths, loops and basic dynamical variables 

Definitions 

Our discussion will start from a space-time of the form X x 1Z, with X being a 
three dimensional surface of constant time, where we have defined the Hamiltonian 
version of our gauge theory. All curves to be handled in the sequel will lay on this 
surface. 

Let us consider a smooth and continuous function r\ : s — > X, where s is a 
parameter in the real interval Si < s < s 2 . An open path on the surface, joining 
the points A (initial) and B (final point) of it, is defined by t)a b (s) = ^ Sl S2 (s) = r), 
with r](si) = A and 77(32) = B. We will use greek letters in the middle of the 
alphabet to denote open paths. The choice to move along the path, either from 
A to B or viceversa, defines its orientation . Given any oriented path t)a b (s), its 
inverse (rj~ 1 )B A (t) is the original path traversed in the inverse sense, with t being 
the corresponding parameter. 

A loop, or a closed path, is an open path where the initial and final points 
coincide, i.e. such that A = B. They will be denoted by greek letters from the 
beginning of the alphabet: a(s), f3(s), for example. Again, the sense of traversing 
the loop defines its orientation. In this way, the inverse of any loop is the same 
loop run over in the opposite sense. Finally, it will prove convenient in the sequel 
to introduce the idea of a multiloop 7 as a collection of loops a±, a 2 , ■ ■ ■ , a p , which 
will be denoted by 7 = a.\ U a 2 U . . . U a p . A multiloop is a direct product of loops 
and does not involve any composition property. 

Composition properties 

It is convenient to define the multiplication of paths, which we denote by o. Two 
open paths i]a b and <Jc D can be multiplied only if either (i) A = D, or (ii) C = B. 
In the first case, 

the composed path, called (a o rf) c B , is the result of going first through the 
path a and subsequently through the path r], in that order. In the second situation 
we obtain the path (77 o (j)a d ■ The multiplication i]a b (ji~ v )b A produces what 



is called a tree, which amounts to start from point A and going back to it, without 
enclosing any area on E. Any two open paths or loops differing by a collection of 
trees, are considered to be equivalent. 

Two oriented loops, a and (3 can be multiplied only when they have a point of 
intersection, say C. Then, the loop a o f3 is obtained starting from the common 
point C, going first through the loop a to end up in C, subsequently going through 
the loop (3, to finally end up again at the intersection point. In each case, the sense 
of travelling is defined by the orientation of the loops to be multiplied. From now 
on we will consider that all loops are based on a fixed point on E, so that they 
always have at least this point in common. This loop space has the structure of a 
semigroup. 

Parallel transport matrix 
To any open path 77 e S, we associate the group element 

U(rj Sl S2 ) = P exp I** ds rfis) Ms), (37) 

J si 

which is the non-abelian generalization of the Aharanov-Bohm phase factor. Here, 
fj z (s) denotes the tangent vector of the path rj and P stands for the path-ordering 
operator defined as 

P exp f 2 dt M(t) = 1 + f 2 dt M(t) + [** dt f* dt x M{t)M(t x ) 

Jsi J Si J Si Jt 

+ r dt r dh f 2 dt 2 ... M(t)M(t 1 )M(t 2 )...:=V(s 1 ,s 2 ), (38) 
Jsi Jt hi 

where M(t) denotes any element of the corresponding Lie algebra, valued along the 
path rj. Let us emphasize that, in general, [M(ti), M(t 2 )} 7^ for arbitrary points 
in the path. The ordering in (38) is defined by si < t < t\ < t 2 < . . . < s 2 . An 
alternative way of writing the generic ordered integral is in the form 

f 2 dt 2 f 2 dh f 1 dt ... M(t)M(ti)M(t 2 ) ■ ■ ■ ■ (39) 

J Sl J Si J Si 

The parallel transport matrix (37) is also known as the integrated connection along 
the corresponding path. As a function of the end points, the group element V(si, s 2 ) 
satisfies the properties 

^fl = V( Sl ,s 2 )M(s 2 ), ^^ = -M( Sl )V( Sl ,s 2 ), (40) 

which can be directly seen from Eq. (39) and Eq. (38) respectively. Also we have 
the composition property 



U( Vsi s >) U{<j S2 s >) = U(( V o a) Sl s *), (41) 

provided the final point of the path 7] coincides with the initial point of the path 
a. In this way, both paths can be composed as indicated in the RHS of Eq. 
(41). The composition appearing in the LHS of Eq.(41) corresponds to the group 
multiplication. Finally, under non-abelian finite gauge transformations, (generated 
by the group element g(x)), Ai — > Ai = gAig^ 1 — gdig^ 1 , the integrated connection 
transforms as 

U(V S1 S2 ) - U(V S1 S2 ) = gisjUiri^g- 1 ^). (42) 



Holonomies and Wilson loops 



The final goal of this formulation is to define strictly gauge invariant objects, 
which will constitute the appropriate variables to formulate the theory. To this 
end, let us consider the holonomy U(a) associated to the loop a, which is just the 
integrated connection (37) around the loop, 



U(a) — P exp f dsa 1, (s)Ai(s). 

J a 



(43) 



This object is still not gauge invariant, but transforms according to (42) with the 
corresponding group elements evaluated at different values of the parameter, which 
nevertheless describe the same point on E: the chosen initial (s = s±) and final point 
(s = S2) of the loop. Group elements must be single valued on E, so that g~ 1 (s 2 ) = 
(g(si))^ 1 . In this way, TrU (a) is indeed a gauge invariant complex number. This is 
the original Wilson loop, whose definition has been generalized to include insertions 
of the gauge- covariant chromoelectric field, as a way to incorporate the variable 
conjugated to the connection in a gauge invariant manner. Then, a tower of Wilson 
loops can be constructed as follows 



TV) 
T\a) 

T ij (a) 



= Tr U(a) = h(a,A), 
= Tr [u(a s ^)E\ Sl )U(a Sl s ») 

= Tr \u(a s ^)E l ( Sl )U(a Sl S2 )E^s 2 )U(a S2 S0 ) 



(44) 



where So parametrizes the origin of the loop a. Each of the numbers in (44) is 
a gauge invariant quantity and we will denote them generically by T° ,T l ,T 2 , . . . 
according to the number of insertions introduced. These objects constitute the fun- 
damental dynamical variables used to formulate the gauge theory in this approach. 

The original symplectic structure (16) of the gauge theory will induce a Poisson 
brackets structure among the Wilson loops of Eq.(44). In a very sketchy form, 
these Poisson brackets will be of the form 



m+n—1 




PB ~ 



(45) 



k=i 



The result in the RHS will involve some recomposition and rerouting of the initial 
loops appearing in the LHS. The details are given in Refs. [15], [17]. 



The dynamical variables defined above , i.e. the Wilson loops, have support 
on the space of closed loops on S, which now replaces the configuration space 
associated to the connection representation. To complete the formulation of the 
gauge theory in the loop space is still necessary to rewrite the Hamiltonian and 
the constraints in the new representation. This can be accomplished through the 
introduction of differential operators which act on loop dependent functions: the 
loop derivative and the connection derivative, for example. The reader is referred 
to Ref. [21] for the details. 

Accordingly, when quantizing the theory, the corresponding operators T n — > T n 
will have also support on the space of loops and they will be constructed as operators 
acting on a Hilbert space which now will be spanned by functionals of loops. An 
heuristical, but very illuminating way of looking at some properties of the loop 
space is via the loop-transform introduced in Refs. [17], 

[18]. This transforms provides the change of basis from the connection represen- 
tation {\A >} to the loop representation {\a >} in the form 



where a is a multiloop with components ai, . . . , a p . The measure [oL4] in (46) is 
not presently known and this fact is what makes this transform only an heuristical 
tool. The wave function of the system |\& > can then be projected either into the 
connection representation < A\ty > or into the loop representation < a\ty >, with 
the two projections being related by the loop transform (46). The change of basis 
(46) will provide also the representation of the operators acting in the loop space, 
starting from their counterparts in the connection representation. 

Nevertheless, there is an alternative way to this construction, based upon the 
formulation of the quantum theory directly in the loop space. To this end, one 
starts from the Poisson brackets relations (45) among the Wilson loops. Following 
the standard procedure, the classical variables are promoted into operators T n and 
their Poisson brackets are translated into commutators. A basis in the loop space 
is constructed by first defining the zero- loop state |0 > such that T n \0 >= 0, n > 
1. An arbitrary multiloop state \a > is subsequently constructed by the action 
T° (cxi)T° (a 2 ) ■ ■ .T°(a p )\0 >:= \a >. Then, using the commutation relations it is 



4.2. The loop space 




(46) 



possible to figure out the action of any operator upon the multiloop states. Ref. 
[16] contains a very nice and explicit application of these ideas. 

Being gauge invariant, the Wilson loop variables automatically satisfy the Gauss 
law constraints of the theory. In this way, the arbitrary multiloop quantum states 
constructed by the action of the operators T° (aii)T ((X2) ■ ■ ■ T°(a p ) on the zero loop 
state, will be automatically annihilated by the quantum Gauss law constraints. 

The equivalence between the connection representation and the loop represen- 
tation of gauge theories defined over a Lie algebra was proved in Ref. [25] for the 
group SU (N). There, it is given a procedure to reconstruct, up to a gauge transfor- 
mation, the gauge field at each point ( written as a complex matrix ) starting from 
the knowledge of all the Wilson loops T°(a) going through that point, together with 
the imposition of the so called Mandelstam identities, which set constraints among 
the Wilson loops. This construction demonstrates that the Mandelstam identities 
are a sufficient set of conditions on the Wilson loop variables which guarantee that 
such variables represent a local gauge theory over a Lie algebra. 

5. OVERCOMPLETENESS OF THE LOOP SPACE 

VARIABLES 

Let us summarize what has been done by saying that the loop space formulation 
of classical gauges theories starts from an infinite set {T n } of numerical, non-local, 
gauge invariant degrees of freedom, defined over loops in the three-dimensional 
surface E. These variables replace those of the standard local connection represen- 
tation, E(x, t), A(x, t), which are not gauge invariant. 

The loop space degrees of freedom are basically traces of group elements, which 
in the sequel will be thought of as traces of matrices of a given dimension. Now, 
it is well known that the Cayley-Hamilton theorem of linear algebra, for example, 
sets relations among the traces of powers of a given matrix, which depend on the 
dimension of the representation. As we will show in the sequel, a similar mechanism 
induces relations among the different Wilson loops T°(a), with the consequence that 
{T n } is an overcomplete set of dynamical variables. These relations among powers 
of traces of products of matrices are generically known as Mandelstam identities 
(MI). Even in the standard case (i.e. gauge theories over a Lie algebra) there is 
lacking a general procedure to obtain the full set of MI for arbitrary restricted 
groups. The situation in the case where the connection is valued on a super Lie 
algebra is still more open and we hope that some of the ideas presented here could 
provide a starting point to deal with these problems. 

At the quantum level, these MI will generate, via the loop transform (46) for 
example, relations among the multiloop states, which will make this set also over- 
complete. 

In this way, any loop space formulation of a gauge theory will require the knowl- 
edge of the corresponding complete set of MI, in order that the true degrees of 
freedom of the system can be identified. The subsequent problem of imposing 



and/or solving the MI is still highly involved and has been carried over, for ex- 
ample, in the following cases: (i) SU(2) pure lattice gauge [26], [27]. (ii) 2 + 1 
dimensional Einstein gravity [28] and (iii) a specific sector of 2 + 1 supergravity 
[29]. 

The observation made above concerning the relation between the MI and the 
Cayley-Hamilton identities will be made more explicit in what follows, by formu- 
lating an algorithm to go from one type of identity to the other, in the generic case 
of matrices over the complex numbers. This procedure will be next extended to 
generate the generic MI for supermatrices, after having also extended the Cayley- 
Hamilton theorem to this case. 



5.1. The Mandelstam identities 

The simplest MI (known as MI of the first class) is just the statement of the 
cyclic property of the trace. For any two complex square matrices M and N, we 
have 

Tr (MN) = Tr (NM) . (47) 

Let us consider the effect of this identity upon the loop space variable T°(a o (3) for 
example. We have T°{a o (3) = Tr U(a o (3) = Tr {U{a)U{(3)) = Tr {U{(3)U{a)), 
where we have used the composition property (41) for the holonomies. The above 
chain of equalities leads to the conclusion 

T°{ao(3) =T°{(3oa). (48) 

When translated into an statement over the quantum loop states, Eq. (47) implies 

\ao(3>=\(3oa>, (49) 

which arises either from the loop transform (46), or from the method of creating 
the quantum loop space vectors starting from the zero loop state |0 >. The above 
constitutes a very simple example of the overcompleteness of both the classical loop 
space variables together with the corresponding quantum states. 

All the remaining MI are known as Mandelstam identities of the second class 
and there is a hierarchy of them, according to the restrictions satisfied by the 
group under consideration. Suppose we are dealing with a group representation in 
terms of n x n matrices M. There will be what we call a generic MI, related to 
the specific dimension n. If we impose further restrictions upon the representation, 
like unit determinant or unitarity, for example, new MI reflecting them have to be 
generated. Most of what we will discuss in these notes has to do with the generic 
MI and their extension to the supersymmetric case. 



The generic Mandelstam identities 



Let us consider the group elements U(a) realized as n x n matrices with no 
further restrictions. The corresponding generic MI are obtained starting from the 
identity [25] 

= e ili2-inin+i e jlj2-jnjn+l ~ ( _ 1) ^ ' ^lPfjl) ' ' ' ^!n+l P 0'n+l)) (^0) 

Pes n+ i 

where ik, jk — 1, • • • n, k — 1, . . . , n, n + 1. Here, the sum is over all permutations 
P of the symmetric group S n+ i and ir(P) denotes the parity of the permutation. 
The zero in the LHS of Eq.(50) comes about because there must be a repetition 
when distributing n objects among (n + 1) places. The second equality is just the 
determinantal expansion of the product of two e-symbols. 

By saturating the relation (50) with n + 1 different matrices one obtains an iden- 
tity among products of traces of products of matrices. The resulting Mandelstam 
identity can be written as [30], [31] 

£ (-lY^WiM,, M 2 , . . . , M n , M n+l ) = 0, (51) 

Perm(l,2,...,n,n+1) 

where 

W(M U M 2 , . . . , M n , M n+1 ) = Tr(M a{1) . . . M a(l) )Tr(M a{t+1) ...)..., (52) 

corresponds to the cycle decomposition (a(l) . . . a(i))(a(i + 1) ...)... of the permu- 
tation P. The generic MI is a function of degree one in each of the n + 1 different 
matrices involved. As an example, let us consider the 2x2 case, which leads to 
the following MI among three matrices 

= Tr(M 1 )Tr(M 2 )Tr(M 3 ) + Tr(M 1 M 3 M 2 ) + Tr(M 1 M 2 M 3 ) 

-Tr(M 1 )Tr(M 2 M 3 ) - rr(M 2 )rr(MiM 3 ) - Tr^M^Tr^M^). (53) 

The restriction implied among the classical loop variables is 

= T°(ai)T (a 2 )T (a 3 ) + T°(«i o a 3 o a 2 ) + T°(ai ° <* 2 ° a 3 ) 

—T°(ai)T°(a 2 o a 3 ) - T°(a 2 )T ( ai o a 3 ) - T G {a 3 )T°{ ai o a 2 ), (54) 

which has the following expression in terms of the quantum loop states 

= \a\ U a 2 U a 3 > +\oti o a 3 o a 2 > +\a\ o a 2 o a 3 > 

\ai U a 2 o a 3 > +\a 2 U cti o a 3 > +\a 3 U oti o a 2 > . (55) 



The restricted Mandelstam identities 



To the author's knowledge, there is no systematic procedure to construct the MI 
that reflect further restrictions on the group elements. In this section we give a 
brief description of the method used in Ref. [31] to implement the MI for the group 
SU(N) realized in terms of NxN matrices Mj in the adjoint representation. In this 
reference, the basic objects to be considered are the numbers Wk{Mi, M 2 , . . . M k ) 
which are defined through the following recurrence relation 

W 1 (M 1 )=Tr(M 1 ), 
(k + l)W k+ i(M u M 2 , . . . M k+1 ) = WiCMfc+OW^Mx, M 2 , . . . M k ) 

-W k (M 1 M k+1 ,M 2 ,...M k ) 
- ... - W k {M u M 2 , . . . M k M k+1 ). (56) 

We have changed the notation with respect to Ref. [31], in order to be consistent 
with our own conventions. In this construction, the generic MI is written as 

W N+1 (M u M 2 ,...M N+1 ) = 0. (57) 

Unitarity is imposed by demanding 

W^M' 1 ) = Wi(M)*, |Wi(M)| < N. (58) 

Finally, the unit determinant condition can be written as Wat(M, M, . . . M) = 1. 
This, in turn, can be reexpressed in terms of N + 1 matrices in the form 

W N (M 1 M N+1 , .U-.U.v.l . . . M N M N+l ) = W N (M U M 2 , . . . M N ). (59) 

This last form is useful to produce a geometric interpretation of the condition of 
unit determinant in the language associated with the loop description of this gauge 
theory. 

5.2. The relation between the Mandelstam and the 
Cayley-Hamilton identities 

In order to motivate the results of this section, let us consider an interesting 
particular case of Eq.(53), by setting \l\ \ l> M, M 3 = X, where X is an 
arbitrary 2x2 matrix. Substituting in Eq. (53) and using the property that 
Tr(CX) = 0, VX, implies C = 0, we obtain the matrix relation 

M 2 - Tr(M)M + det(M)l 2 = 0, (60) 

where X 2 is the 2x2 identity matrix. In the above we have used the fact that 
det(M) = \ (Tr(M) 2 - Tr(M 2 )) for 2 x 2 matrices. Next, we recognize that 
Eq.(60) is nothing but the statement of the Cayley-Hamilton theorem for the 2x2 
matrix M. As we will show in the sequel, this construction can be generalized to 
n x n matrices and can also be used in the reversed sense, that is to say, starting 
from the identities arising from the Cayley-Hamilton theorem, an algorithm to 
obtain the generic MI in the general case can be constructed. 



The Cayley-Hamilton theorem in classical linear algebra. Let us pro- 
vide a brief review of this important theorem. For an arbitrary n x n complex 
matrix M, the corresponding characteristic polynomial P(x) is defined as 



P(x) = det (xl n - M) =x n + a lX n ~ l + ...a n . (61) 

The Newton equations provide a recursive method to calculate the coefficients <2j, 
in terms of traces of powers of M [32] 

1 ' 

o-i+i = ~ t— : ij+i-fc o,k, (62) 
where t/~ = Tr(M k ). An explicit solution of the above recursion relations is 

ao+Ql+ ... +as=i+ i ("o + «i + • • • + a s )(ai + . . . + a s ) . . . {a s ) 

where the summation is over all the unordered distinct partitions (ao, . . . , a s ) of 
i + with s + 1 being the total number of terms in each partition. 

Extending the complex variable x to a matrix-valued variable, the Cayley- 
Hamilton theorem states that each matrix annihilates its characteristic polynomial 
[33], i.e. 

P(M) = M n + a 1 M n ~ 1 + ...a n l n = 0. (64) 

Let us emphasize that (64) is a matrix identity involving n 2 numerical relations. 
We will refer to the relations of the type (64) as the Cayley-Hamilton identities 
(CHI). 

Now we describe the proposed algorithm to obtain the MI starting from the 
CHI [34]. The basic idea is to produce a sequence of CHI that leads to the MI, 
incorporating n + 1 different matrices and recalling that the final MI must be an 
homogeneous function of degree one in each matrix. The procedure is as follows: 
starting from the identities P Ml { M i) = °, Pm 2 {M 2 ) = 0, P Mi+ m 2 { m i + M 2) = 
we construct 

T 2 (Mi, M 2 ) := P Mi+ m 2 (Mi + M 2 ) - P Ml (Mi) - Pm 2 {M 2 ) = 0. (65) 

We are denoting by Pm{x) the characteristic polynomial corresponding to the ma- 
trix M and also we have set Pm{M) := 7\(M). Every term in (65) is an homoge- 
neous function of Mi and M 2 and the subtractions are designed in such a way that 
both Mi and M 2 appear at least once in every term of T 2 (Mi,M 2 ). In this way 
T 2 (Mi = 0,M 2 ) and T 2 (M 1: M 2 = 0) are identically zero. Moreover, we consider 
that T 2 (Mi, M 2 ) can be fully expanded using the distributive property of both the 
trace and the matrix product with respect to matrix addition. The next step is to 
construct 



T 3 {M 1 ,M 2 ,M 3 ) : 

red: 

(Mi + M 2 + M 3 ) - T 2 (Mi, M 2 ) - T 2 (M ± , M 3 ) 
-T 2 (M 2 , M 3 ) - T 1 (M 1 ) - T 1 (M 2 ) - Ti(M 3 ) = 0. (66) 

Again, T 3 is identically zero whenever any of the Mj is set equal to zero. We have 
introduced the subscript \ re( i to indicate an identity which has been reduced in 
such a way that every matrix involved is present at least once in every term, after 
the identity is fully expanded. In other words, the expression Pm 1 +m 2 +m 3 (Mi + 
M 2 + M 3 )| re d in (66) can be directly constructed by expanding the corresponding 
characteristic polynomial and discarding all terms in which any one of the three 
matrices is missing. 

Extending this idea, we construct reduced identities of always increasing order, 
where we subtract all the lower order identities at our disposal. This procedure 
leads to 

T k {M 1 ,...,M k ) : = PM 1 +...+M k (Mi + . . . + M k ) 

-Y,T i (M 91 ,...,M ai ), (67) 

i<k 

where the sum is carried over all subsets {si, . . . , Sj} of {1, . . . , k}. The fact that 
the characteristic polynomial is of order n guarantees that T k (Mi, . . . , M k ) are 
identically zero for k > n + 1. The expression T n (Mi, . . . , M n ) is homogeneous of 
degree one in each of the n matrices involved. 

The generic MI is obtained from Tr (T n (Mi, . . . , M n )M n+ i) as shown with all 
detail in Ref. [34]. The converse statement, which amounts to recovering the CHI 
(64) starting from the MI (51), is also proved in general in this work. 

Before closing this section let us observe that the immediate application of the 
above method to the construction of the MI in the case where further restrictions 
upon the group elements are required, will work only if such restrictions are pre- 
served by matrix addition [35]. In particular, the restriction to matrices with unit 
determinant cannot be directly implemented in this way. It remains still an open 
problem to work out a general procedure to construct the MI corresponding to 
arbitrarily restricted groups. 

Nevertheless, there are two interesting examples where a restricted MI is directly 
obtained from the characteristic polynomial. These correspond to the cases of the 
groups SL(2,$l) and SU(2), which are relevant in the study of 2 + 1 de Sitter 
gravity, and 3 + 1 Einstein gravity in the Ashtekar formulation, respectively. Let 
us consider two of such 2x2 matrices, Mi and M 2 , with unit determinant. The 
CHI for the first matrix can be rewritten as 

Mi - Tr(Mi)J 2 + Mf 1 = 0, (68) 

after the Eq. (60), with det(Mi) = 1, has been multiplied by Mf 1 . Multiplying 
further Eq. (68) by M 2 and taking the trace, we end up with the restricted MI 



Tr(MiM 2 _1 ) + Tr(M 1 M 2 ) = Tr(M 1 )Tr(M 2 ) 



(69) 



which involves only two matrices, instead of three as in the generic case presented 
in Eq. (53). 

In the case of the group SL(2,$t) it is also possible to exhibit an ex- 
ample of the drastic reduction of the independent loop space variables pro- 
duced by the ML Let us consider the infinite set of Wilson loops of the form 
Tr ( Mf 1 Mf 1 Mf 2 Mf 2 . . . Mf"Mf™ . . .), for any integer Pi , ft . Using the MI (69) one 
can show that any of such Wilson loops can be expressed as a function of three 
traces only: Tr(M 1 ), Tr(M 2 ) and Tr(M l M 2 ) [36]. Thus, in the corresponding 
sector of the theory we will have only three independent degrees of freedom. A 
simple example of such reduction is to consider Tr(MfM 2 ) for example. Here we 
apply the relation (69) with Mi -> M u M 2 -> M X M 2 obtaining 



Our main motivation for the work described in the previous section has been 
the possibility of extending these ideas to the loop space formulation of pure gauge 
theories having their connection valued on a super Lie algebra, as it is the case of 
supergravity in the Ashtekar variables formalism [8], for example. In this situation, 
the group elements are described by generic supermatrices as defined in Eq. (28). 
The knowledge of the MI for supermatrices will be also relevant to the loop space 
formulation of fully supersymmetric gauge theories. 

It is well known that Grassmann numbers can be realized in terms of complex nu- 
merical matrices. From this point of view one could think that the Cayley-Hamilton 
theorem for supermatrices would be just a trivial extension of the standard case. 
Nevertheless, this is not the case for at least two reasons: (i) after a realization of 
the Grassmann numbers in terms of numerical matrices it would not be possible 
to recover a result in terms of the original Grassmann numbers . (ii) the matrix 
realization will effectively augment the size of the resulting supermatrix , which 
now would be completely numerical, in such a way that the resulting standard 
characteristic polynomial would be also of higher degree. 

Let us give a precise meaning to the above observations using the simplest case 
of a (1 + 1) x (1 + 1) supermatrix 



Tr[M\M 2 ) = Tr(Mi)Tr(MiM 2 ) - Tr(M 2 ). 



(70) 



6. THE CAYLEY-HAMILTON THEOREM FOR 
SUPERMATRICES 




(71) 



where p, q are even Grassmann numbers while 77, 9 are odd Grassmann numbers 
such that rf = = 6 2 , rj9 = —Or]. The minimum size for the 



realization of these Grassmann numbers corresponds to 4 x 4 gamma-matrices, 
like i] = 7o + 7i, = 72 + 273 in the signature (+, — , — , — ), for example. In this way, 
M 2 will be realized as an 8 x 8 numerical matrix and, consequently, the resulting 
characteristic polynomial will be a numerical polynomial of degree 8. Clearly, it will 
not be possible to rewrite the numerical coefficients of the polynomial in terms of 
the original Grassmann numbers. Nevertheless, it is possible to find a characteristic 
polynomial of degree 2 for this case, which is given by 



We can verify, by direct substitution, that the above polynomial satisfies P 2 (M 2 ) = 
0. In this simple case, the polynomial (72) can be constructed just by solving this 
condition. When (q — p) 2 7^ 0, the polynomial (72) can be redefined in the monic 
form of Eq.(61) with the choices a\ = q + p — -^~^]0 and a 2 = qp — ^z^V®- Thus, a x 
generalizes Tr(M 2 ) = p + q, while a 2 generalizes det(M 2 ) = qp — r\0 corresponding 
to the case where all entries in M 2 would be complex numbers. Here we also see 
that the standard determinant is not well defined for supermatrices: we would 
have to make a choice among all possibilities qp — (Ar/9 + (1 — A)0rf), for arbitrary 
A. In this way, the heuristical result (72) certainly motivates the search for a 
general procedure to construct such null polynomials, leading to the construction 
of the Cayley-Hamilton theorem for supermatrices. This is an interesting problem 
in its own, besides the possible applications to the loop space formulation of gauge 
theories involving supersymmetry. 

6.1. The characteristic and null polynomials for 

supermatrices 

The first step in this direction is to provide a definition of the characteristic 
polynomial, which extends the standard one given in Eq. (61). This problem is 
certainly related to the eigenvalue problem of a supermatrix, which is discussed 
in Ref. [37]. The eigenvalues of an (m + n) x (m + n) supermatrix M are even 
Grassmann numbers which can be of two types: (i) first-class eigenvalues, Xi, i = 
1, . . . , m , whose eigenvectors are of the form [Ei, . . . , E m , Oi, . . . , O n ] T , where Ek 
denotes even Grassmann numbers, while 0/ labels odd Grassmann numbers, (ii) 
second-class eigenvalues, A a , a = 1, . . . ,n with corresponding eigenvectors of the 
form [Oi, . . . , O m , Ei, . . . , E n ] T . The characteristic polynomial will certainly read 



in terms of the eigenvalues of the supermatrix . When written in the form of Eq. 
(61), the explicit expressions for the coefficients will be given again by Eq. (63) 
with t k — Xi + . . . + A^j + Xi + . . . + A£ . A problem now arises if we want to rewrite 
such coefficients in terms of the supermatrix itself : the appropriate invariants for 
this case are not traces of powers of the supermatrix, but supertraces of powers of 



-fM^) = (<? — p)x 2 — (q 2 — p 2 — 2r)0)x + ((q — p)qp — (q + p)rj0). 



(72) 



P(x) = (x - Xi) . . . (x - X m )(x - Xi) . . . (x - A n ), 



(73) 



the supermatrix instead. The latter are given by Str(M k ) = \\ + . . . + A^ — A^ — 
. . .-A*. These expressions break the permutation symmetry among all eigenvalues, 
leaving only a permutation symmetry among the first-class eigenvalues together 
with an independent permutation symmetry among the second-class eigenvalues. 
The coefficients are nevertheless symmetric in the whole set of eigenvalues and 
this is precisely why it is a more involved task to rewrite them in terms of the 
corresponding supertraces. 

It would seem that a reasonable starting point for the construction of null poly- 
nomials in the case of supermatrices is the natural extension of the definition (61) 
to P(x) — > Sdet(xT — M). Nevertheless, this function is not a polynomial but a 
ratio of polynomials, as can be seen from its expression in terms of the eigenvalues 

^*-*-%--t:%~-KY (74) 

As a preliminary step towards the definition of the null polynomial associated to 
the supermatrix (28), let us introduce the standard null polynomials corresponding 
to the even block-matrices A and D 

a(x) = det(xX m — A), d(x) = det(xl n — D), (75) 

respectively. Next, let us consider the general expression for 

h(x) := Sdet(xl m+n - M), (76) 

which will be called the characteristic function in the sequel. From the two alterna- 
tive expressions to calculate the superdeterminant, given in Eq. (32), we are able 
to write the characteristic function as 

where the basic polynomials F, G, F and G are given by 

F(x) = det{d{x){xl -A)- Badj{xl - D)C), G{x) = (d(x)) m+ \ 

F(x) = (a(x)) n+1 , G(x) = det(a(x)(xl — D) — Cadj(xl - A)B). (78) 

The above expressions are readily obtained from Eqs. (32), using the relation 
(xl — F) _1 = [det(xl — F)]~ 1 adj(xl — F) valid for any even matrix F. Notice that 
F is expressed in terms of the determinant of a m x m even matrix, while G(x) is 
the determinant of a n x n even matrix. 

In order to motivate the basic idea of our definition for the characteristic polyno- 
mial of a supermatrix, let us consider the simple case of a block-diagonal superma- 
trix M (i.e. B — 0, C — 0). Here h(x) = a(x)/d(x) and clearly the characteristic 
polynomial is P(x) = a(x)d(x), which is the product of the numerator and the 
denominator of the corresponding superdeterminant. In fact we have 



because a(A) = 0, d(D) = 0. An analogous statement is obtained for a supermatrix 
written in terms of its eigenvalues. In the general case, where h(x) is given by 
Eq.(77), the numerator of the superdeterminant is F (F) while the denominator 
is G (G). These observations lead to the following definition of the characteristic 
polynomial of a supermatrix 

V(x) := F(x)G(x) = F(x)G(x), (80) 

in terms of the basic polynomials F, G, F and G, given in Eqs.(78). For notational 
simplicity we will not necessarily write explicitly the re-dependence on many of the 
polynomials considered in the sequel. 

Let us consider again the block-diagonal case, this time when a(x) and d(x) have 
a common factor f(x), i.e. 

a(x) = f(x)a 1 (x), d(x) = f(x)d 1 (x). (81) 

In this example, a null polynomial is given by P(x) = f(x)ai(x)di(x), which is 
a polynomial of degree lower than the product a(x)d(x). Motivated by this fact, 
together with the work of Ref. [38], we realize that there are some cases in which 
we can construct null polynomials of lower degree than V(x), according to the 
factorization properties of the basic polynomials F, G, F, G. 

At this point it is important to observe that we do not have a unique factorization 
theorem for polynomials defined over a Grassmann algebra. This can be seen, for 
example, from the identity x 2 = (x + za)(x — za), where a is an even Grassmann 
with a 2 = and z is an arbitrary complex number. 

The construction of the null polynomials of lower degree starts from finding the 
divisors of the pairs F, G, (F, G) which we denote by R, (S) respectively. This 
means that one is able to write 

F = Rf, G = Rg, 

F = Sf, G = Sg, (82) 

where all polynomials are monic and also /, g, f, g are of lower degree that their 
parents F, G, F, G, by construction. They must satisfy 

f/g = f/g, (83) 

because of Eq. (77). The expressions in (82) might not be unique. 

Let us emphasize that in the case of polynomials over the complex numbers, when 
R and S are of maximum degree, Eq. (83) would imply at most / = \f, g — Xg 
with A being a constant. Since we are considering polynomials over a Grassmann 
algebra, this is not necessarily true as can be seen again in the above mentioned 
identity xj[x — za) — (x + za)/x, which we have rewritten in a convenient way. 



The above discussion leads to the following definition: given an arbitrary (m + 
n) x (m + n) supermatrix M, with a characteristic function h(x) such that the 
polynomials F, G have a common factor R and the polynomials F, G have a common 
factor S, satisfying Eqs. (82) and (83), then a null polynomial of M is given by 

P(x) := f(x)g(x) = f(x)~g(x). (84) 

The polynomial (84) is clearly of lower degree than V(x), which is just a particular 
case of the null polynomials (84) when R — S — l. We will concentrate mostly on 
(84) in the sequel. 

6.2. Proof of the Cay ley-Hamilton theorem for 

supermatrices 

In this section we show that the polynomial defined in Eq.(84) does in fact 
annihilates the supermatrix M. To this end, we first extend a lemma often used 
to prove the Cayley-Hamilton theorem for ordinary matrices [33]. We briefly recall 
such lemma and emphasize that it is independent of the matrix considered being a 
standard matrix or a supermatrix. It goes as follows: let M, (xl — M) and N(x) 
be (m + n) x (m + n) supermatrices, where M is independent of x. Let N(x) be 
a polynomial supermatrix of degree (p — 1) in x, i.e. N(x) = iV a; p ~ 1 + Nix p ~ 2 + 
.. + N p _ix°, (where each N k , k — 0, • • • ,p — 1, is a (m + n) x (m + n) supermatrix 
independent of x), such that 

(xl m+n - M)N(x) = P(x)l m+n , (85) 

where P(x) = a x p + aix^ 1 + • • • + a n x° is a numerical polynomial of degree p. 
Then, one can prove that P(M) := a M p + cqikP" 1 H — • + a n T m+n = 0. The proof 
follows by comparing the independent powers of x in Eq. (85) and then explicitly 
computing P(M) [33]. 

In the standard case, the polynomial matrix N(x) is just given by N(x) = 
adj(xl - M) = det(xl - M)(xl - M)" 1 , and P(x) = det(xl - M). In the case 
of a supermatrix we do not have an obvious generalization either of the matrix 
adj(xl — M) or of det(xl — M). Nevertheless, following the analogy as close as 
possible we define 

N(x) := P(x)(xl - M)-\ (86) 

where P(x) is the polynomial introduced in Eq.(84) of the previous section. 

The challenge now is to prove that N(x), which trivially satisfies the Eq. (85), 
is indeed a polynomial matrix. In this way we would have proved that P(M) = 
according to the property stated after Eq.(85). 

To begin with, we show that the blocks corresponding to the inverse supermatrix 
(xX m+n — M)~ l can be written in a compact form as 



1 flip 1 Q/~< 

(xl m+n - M)~ l = (xX m+n - M)" 1 = -_ (87) 

1 OF 1 dG 

(xI m+ „ - M)~] = Jm -, (xl m ,„ - M)~l = (88) 

where Aij, Bj a ,C a jand D a/3 are the entries of the supermatrix M defined in Eq. 
(28) and F, G, are the polynomials given in Eqs. (78). The derivative with respect 
to an odd Grassmann number is taken to be a left derivative, defined such that 
5F = 5Bj a -^—. The proof of the above equations begins with the calculation of 
(xX m+n — M)" 1 in block form, with the results 

(xX m+n - M)- 1 = ((xX m -A)- B(xX n - D)- l C)-\ 
(xX m+n - M)- 1 = -(xX m - A)- l B{{xX n — D) — C(xX m - A)~ 1 B)~ l , 
(xX m+n - M)zi = -(xX n - D)- l C{{xX m -A)- B(xX n - D^C)' 1 , 
(xX m+n - M)- 1 = ((xX n — D) — C(xX m - A)- l B)-\ (89) 

Here, the subindices 11, 12, 21 and 22 denote the corresponding mxm,mxn,nxm, 
and n x n blocks. The above block form in Eq. (89) has the same structure as 
in the classical case. Let us concentrate now in the 11 block. Rewriting all the 
inverse matrices of the first Eq.(89) in terms of their adjoints, together with the 
corresponding determinants, we obtain 

(xX m+n - M)ii = tadj((xX m - A)d - Badj(xX n - D)C). (90) 
t 

On the other hand, using the basic property 

5detQ = Tr(adjQ5Q), (91) 

valid for any even matrix Q, we can calculate the change of F with respect to Aij, 
keeping constant all other entries, obtaining 

5F = -d [adj((xX m - A)d - Badj(xX n - D)C] ij 5A jh (92) 

which can be written as 

dF 



-d [adj((xX m - A)d - Badj(xX n - D)C% . (93) 
OAji 

The comparison of Eq. (93) with Eq. (90) completes the proof of the first relation 
in Eq. (87). The proof for the remaining Eqs. (87-88) is carried along similar steps. 

We observe that the conditions for the existence of (xX m+n — M)~ l are the same 
as those for the existence of Sdet(xX m+n — M) which read det(xX m — A) ^ and 
det(xX n —D) ^ 0. Since x is a generic even Grassmann variable, we will assume that 
this is always the case. In this way, the term ((xX m — A) — B[xX n — D)~ 1 C)~ 1 , for 



example, can always be calculated as (X m — [xl m — A) 1 B(xI n — D l )C) 1 (xl m — 
A)^ 1 . The factor on the left can be thought as a series expansion of the form 

1/(1 - z ) = 1 + z + z 2 H , with z = (xl m - A)- 1 B(xl n - D)- l C. Moreover, the 

series will stop at some power because z is a nilpotent matrix. 

Now we are in position to show the principal result of this section which is 
that N(x) := P(x)(xl m+n - M)~\ with P(x) given in Eq.(84), is a polynomial 
supermatrix. Let us consider the block-element 11 of N(x) to begin with. According 
to Eqs. (87-88) together with Eq. (82), this block can be written as 

AT = -n 1 . (94) 

The first term of the RHS in Eq. (94) is clearly of polynomial character. In order 
to see that the second term is also polynomial, we make use of the property 

dlnG g dlnR ^ ding 
dAji dAji dAji 

which follows from the factorization G = Rg, together with the fact that G is just 
a function of D a/ 3, according to the first Eq. (78). In this way, and using also the 
Eq.(83), we obtain 

N " = f M-- g M-- (96) 

which leads to the conclusion that the block-matrix iVV, is indeed polynomial. The 
proof for N ai runs along the same lines, except that now the derivatives are taken 
with respect to B ia and that we have to use |^ = 0, instead of Eq. (95). The 
remaining terms N ia and N a p can be dealt with in analogous manner, by considering 
the derivatives of G = Sg with respect to C a i and Dp a , and by replacing the 
condition (95) by |^ = and = respectively. The results are again of 

the form (96), the only difference been the variables with respect to which the 
derivatives are taken. 

Finally, we can state the following extension of the Cayley-Hamilton theorem to 
the case of supermatrices [39], [40], [41]. 

Theorem (Extended Cayley-Hamilton Theorem) Let M and (xl — M) be (m + 
n) x (m + n) supermatrices, with x being a generic even Grassmann variable. Let 
also Sdet(xX m+n — M) = F/G = F/G, where the polynomials F,G,F and G are 
given in Eqs. (78). Then, for any common factor R such that F = Rf, G = Rg 
and S such that F = Sf, G = Sg, where f/g = f/g, the polynomial P(x) = 
f(x)g(x) = f(x)g(x) annihilates M, i.e. P(M) = 0. 



6.3. Examples of null polynomials for supermatrices 



Here we present two simple examples of null polynomials, constructed according 
to the procedure stated in the last section 



The case of (1 + 1) x (1 + 1) supermatrices. Let us consider the super- 
matrix (71) with d 7^ q. Here the bar denotes the complex component of an even 
Grassmann number, called the body in the literature. From Eqs. (78) we obtain 
the following basic polynomials 

F = (x-q)(x-p)-ri6, G — (x - q) 2 , 

F={x-p)\ G=(x-q)(x-p)+r,9, (97) 
The above functions can be rewritten as follows 

W f M 719 \( ^ \ 

t = [ x — p H \ [ x — q 



q — p ) \ q — p 



q-pj \ q — p / 

r, ( ¥ \ ( V 9 s 

t = [ x — p H \ [ x — p 



q — p J \ " q — p J 

G=(x-q+^)(x-p-^-). (98) 
V Q-P/ \ Q-PJ 

Thus, the factorization properties of Eq.(82) are realized with 

R=L- q -JL\, s=(x- P - 719 ^ 



q — p J \ q — p ) 

f = f=(x-p + —), 9 = 9=(x-q + —)- (99) 

In this way, using the definition (84), we recover the null polynomial of minimum 
degree given in Eq.(72). 

The case of Osp(l\2;C) supermatrices. Another simple example corre- 
sponds to the case of supermatrices belonging to the supergroup Osp(l\2; C), which 
are relevant in the description of de Sitter supergravity in 2 + 1 dimensions [29]. 
We consider a (2 + 1) x (2 + 1) realization of this supergroup defined by the set of 
all supermatrices M which leave invariant the supersymplectic form H 



M r HM = H, H = -1 , (100) 




where T denotes the supertransposed. The supermatrices in Eq.(lOO) can be 
parametrized in the following way 

M= (i i), { = (+Y (WD 



where A is a (2 x 2) even matrix; Xi, x 2 are arbitrary odd Grassmann numbers and 
the superindex T denotes standard transposition. The condition (100) translates 
into the following relations among the elements of Eq. (101) 

X T = £ T EA, a = l + x 1 x 2 , det(A) = 1 - x x x 2 , (102) 

where E denotes the 2x2 antisymmetric block of H in Eq.(lOO). Using the algo- 
rithm of the previous section, we conclude that the unique irreducible expression 
for the characteristic function is 

h(x) = + (103) 

which means / = f , g — g. The null polynomial of minimum degree is then [29] 

P{x) = fg = x 3 - (2 + Str(M))(x 2 - x) - 1. (104) 



6.4. Two examples of Mandelstam identities for 

supermatrices 

As a first step towards the search of an algorithm to produce the MI for super- 
matrices, one may try to directly extend the procedure developed in section 5.3 for 
the case of ordinary matrices. The starting point now will be the null polynomials 
constructed in section 6.1., which can be rewritten in terms of a finite number of 
supertraces. In relation to this, we still suffer a main drawback which is the lack 
of knowledge of a recurrence that would allow to obtain a closed expression for the 
coefficients of the null polynomial in terms of supertraces, in a manner similar to 
the standard case. Another new feature is that the null polynomials are not monic 
any more, in such a way that a would now be a function of the supertraces. This 
property will effectively raise the degree of homogeneity of the supermatrices in 
the corresponding CHI. Nevertheless, such identities will be homogeneous of some 
degree, say t for example. This will allow us to make the following definition of the 
corresponding generic MI 

Str (P Ml +M 2 +...+M t (Mi + M 2 + . . . + M t )\ red M t+1 ) = 0. (105) 

The case of (1 + 1) x (1 + 1) supermatrices. In this situation, the null 
polynomial (72) can be rewritten as 

Str(M)M 2 - (Str(M 2 ))M + ^ (Str(M z ) - Str(M) 3 ) J 2+1 = 0, (106) 

in terms of supertraces, with t = 3. Using the definition (105) we obtain 



Str(A) (Str(BCD) + Str{CBD)) + Str(B) (Str(ACD) + Str{CAD)) 
Str{C) (Str(ABD) + Str(BAD)) + Str(D) (Str(ABC) + Str(BCA)) 
-2Str(AB)Str(CD) - 2Str(BC)Str(AD) - 2Str(AC)Str(BD) 
-2Str(A)Str(B)Str(C)Str(D) = 0, (107) 

which corresponds to a symmetric MI of order four. We have verified this identity 
using Mathematica. 

The case ofOsp(l\2;C) supermatrices. This is an example of a restricted 
MI, which is the supersymmetric analogue of the identity (69), valid for the group 
577(2). Starting from the null polynomial (104) for the supermatrix Mi, multiplying 
this equation by M 1 ~ 1 M 2 and taking the supertrace we are left with 

Str(M 2 Ml) - (2 + StriM^iStr^Mi) - Str(M 2 )) - Str(M 2 M~ 1 ) = 0. (108) 

The above identity has been useful in the identification of the true degrees of 
freedom on one sector of 2 + 1 super de Sitter gravity. For one genus of the 
generic spatial surface, the most general Wilson loop variables are the infinite set 
of supertraces: Str[M{ x Mf M{ 2 Mf . . . Mf n M 2 n . . .), for any integer p u q t . Using 
the MI (108), it is possible to reduce the above infinite set of observables to only 
five complex quantities, which are Str(Mi), Str(M 2 ), Str(MiM 2 ), StrlMiM^) and 
Str(M 1 M 2 MlMl) [29], in complete analogy with 2 + 1 de Sitter gravity [36]. 

7. SUMMARY AND OPEN PROBLEMS 

We have presented a rather sketchy review of the loop space formulation of gauge 
theories, which does not make full justice to all the numerous achievements and 
applications that this method has produced so far. We have tried to incorporate 
a reasonable, but certainly not exhaustive, list of references which remedy this 
situation, offering the reader a detailed version in each situation. 

Our major emphasis has been in the discussion of the formulation of the Mandel- 
stam identities, which appear as unavoidable constraints either among the Wilson 
loop variables or among the quantum loop states, that constitute the natural de- 
grees of freedom of the method in the classical and quantum situation, respectively. 
In the case of pure gauge theories over a Lie algebra, we have shown the equivalence 
between the generic Mandelstam identities, for a given dimension of the representa- 
tion, and the identities arising from the application of the Cayley-Hamilton theorem 
to the matrices of such representation. We have provided an algorithm to go in 
either direction. The main thrust of this development has been it extension to 
provide a general procedure to construct the generic Mandelstam identities in the 
case of pure gauge theories defined over a super Lie algebra, as a first step to deal 
with fully supersymmetric theories. A previous step in this direction has been the 
formulation and proof of the Cayley-Hamilton theorem in the case of supermatrices. 
The extension of this theorem proceeded in two steps: (i) the identification and def- 
inition of a characteristic polynomial for supermatrices and (ii) the proof that each 



supermatrix annihilates the polynomial previously defined. Furthermore, starting 
from the characteristic function (76), we have described a systematic method for 
constructing null polynomials for supermatrices. The resulting Cayley-Hamilton 
identities can be subsequently used to derive the corresponding Mandelstam iden- 
tities. Two simple examples were presented. The construction of the Mandelstam 
identities for arbitrary restricted groups or supergroups remains still an open prob- 
lem. 

Another very interesting question in the loop space formulation of gauge theories 
is the inclusion of fermions, which has been the subject of recent investigations. 
There are at least three ways of approaching this problem: (i) one is to consider 
the fermions as standard matter to be coupled to the gauge theory, (ii) other 
possibility is to consider them as pieces of a superconnection, as has been done 
in these work and (iii) the final possibility is to introduce the fermions in a fully 
super symmetric theory, as partners of integer spin fields. In the first case, several 
matter fields like electrons [42] and quarks [43] have been taken onto account in 
the loop space picture. Also, the introduction of fermions in Einstein gravity has 
been considered in Ref. [44]. The basic idea in these works is to define additional 
gauge invariant variables, besides the Wilson loops, represented by open paths 
which start and end up in the fermions. As mentioned previously in the text, the 
possibility (ii) has been already considered in the case of 2 + 1 super Chern-Simons 
theories and also in the recent discussion of 3 + 1 canonical supergravity in the 
loop space approach. From a general point of view, there is still lacking a proof of 
the equivalence between the proposed loop space representation of these theories 
and the standard connection-matter formulation of them. This will require, among 
many other developments, the construction of generalized Mandelstam identities 
including the open path variables together with the generalization of Giles theorem 
to the superconnection case. Finally, the alternative (iii) has recently being explored 
in Ref. [45], where the fully supersymmetric Wilson loop has been constructed in 
terms of chiral superfields and supercurrents in superspace. 
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